Vibrational mechanics in an optical lattice: controlling transport via potential 

renormalization — Supplemental Material 



A. Wickenbrock, P.C. Holz, N.A. Abdul Wahab, P. Phoonthong, D. Cubero, and F. Renzoni 



o 

(N 
O 

Q 



in 
c3 



c 

o 
o 



(N 



X 



I. MODEL AND DEFINITIONS 

As a model for our experiment, we consider the sim- 
plest configuration in which Sisyphus cooling has been 
shown to take place, i.e. the case of a J g = 1/2 —> J e = 
3/2 atomic transition in a ID optical lattice generated 
by two counterpropagating laser fields with orthogonal 
linear polarizations. This is the so-called lin _L lin config- 
uration [1] . In Rcf. [2| , the following generalized Fokker- 
Planck equation was found in the semiclassical limit for 
the probability density P±(z,p, t) of each atom that is in 
the ground state sublevel |±) = \ J g = 1/2, M g = ±1/2) 
at the position z with momentum p: 



dt m dz ± ^ ^ dp ~^ dp 



P± = 
l±{z)P ± + lT {z)P T 



dp 



— [D ± (z)P ± +L ± (z)P T ], (1) 



where m is the atomic mass and U±(z) = dU±(z) / dz 
U±(z) = T ^-[-2±cos{2kz)] 



(2) 



is the optical bipotential created by the laser fields, with 
k the laser field wave vector; F(t) is a time-dependent 
driving force that can be generated by phase modulating 
one the lattice beams 0]; 



r' 

7±0) = y[l±cos(2fcz)] 



(3) 



is the transition rate between the ground state sublevels, 
with r' the photon scattering rate; 



D ± (z) 



7h 2 k 2 T' 



90 



[5 ± cos(2fcz) 



(4) 



is a noise strength coefficient describing the random mo- 
mentum jumps that result from the interaction with the 
photons without transition between ground state sub- 
levels; and 



L ± {z) 



n 2 k 2 r 



90 



[6 T cos(2&z)] 



(5) 



is related to random momentum jumps that appear in 
fluoresecence cycles when the atom undergoes a transi- 
tion between the atomic sublevels. The normalization 
condition is given by 



J dz J dp[P-(z,p,t)+P+(z,p,t)} = l. 



(6) 



We consider two different types of time-dependent 
driving forces: a bi-harmonic drive of the form 

F d {t) = A x cos(cjf) + A 2 cos(2wi + </>), (7) 

and a high-frequency (HF) drive of the form 

F HF (t) = A HF sin(w H Fi + 4>o)- (8) 

The bi-harmonic drive is used to probe the potential 
amplitude. The HF drive determines the potential renor- 
malization, as discussed in next Section. 



II. POTENTIAL RENORMALIZATION BY 
HIGH-FREQUENCY DRIVING 

We now study the effect on the cold atom system of a 
high-frequency signal Fhf, of the form of Eq. [8] with <j)Q 
an arbitrary initial phase. A low-frequency bi-harmonic 
drive Fd(t) is also included in the analysis, to model the 
experiments in which the potential is probed by using 
this type of driving. 

We are interested in situations in which the frequency 
whf is much larger than uj and any other characteristic 
frequency in the system. Formally, this can be achieved 
by taking the asymptotic limit whf — > oo. In this limit, it 
is also necessary that ^4hf — > oo if the HF signal is to have 
any effect. Due to this strong driving, the momentum 
changes very rapidly, since its time-derivative is of order 
Ahf- Integrating this dominant term in time, we find a 
rapidly changing contribution to the position z(t) that 
goes as — r sin(o;HFi + ^o)j where 



A 



HF 



(9) 



^HF 



Formally, we will consider the asymptotic limit 
whf,^4hf oo while keeping r fixed. By extracting 
the fast dependence from z(t), 



z(t) = z(t) + r sin(w H F^ + <M, 



(10) 



it is expected that z(t) changes on a much slower time- 
scale than that of the HF signal. The density probabili- 
ties for the new variable are then given by 



P±(z,P,t) = P±[z - rsin(w H F^ + <fo), 

p - ™iwhf cos(whf£ + 4>o),t], 



(11) 



where p = mdz/dt. The corresponding generalized 



2 



Fokkcr-Planck equations are 



d p d 



U> ± (z,t)^ + F d (t)?- 
at m oz op op 



dp 2 



-l±{z,t)P ± +%(z,t)P T 
D±(z,t)P ± + L±(z,t)PJ , (12) 



where U' ± (z,t) = U'±[z — r sin(o;HF* + 4>o)], 7±(-M) = 
~/±[z — rsin(a;HFi + <A))L and similarly for D±(z,t) and 
L±{z, t). These coefficients depend on time only through 
the HF signal. On the other hand, both P± and Fj vary 
with time on a much longer timescale. Therefore, we 
could remove the time dependence from the above coef- 
ficients by integrating over a time interval that includes 
many HF periods but in which P± and Fd do not appre- 
ciably change. Equivalently, we can eliminate that fast 
dependence by noting that P± should be independent of 
the HF phase <j)o- Integrating Eq. (fT2"j) over the phase 
</>o, we finally find a generalized Fokker-Planck equation 



analogous to (fT2|) but with the following coefficients: 

,-2k 



U±(z) 



i>oU±[z-r sinfwHFi + <Po)] 



1 

2^ 

= y[~ 2± Jo(2Ax)cos(2fc£)], 

r' 

l±(z) = g-[l ± Jo(2/cr) cos(2fcz)], 



D ± (z) 
L ± (z) 



in 2 k 2 v , 
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where 



(13) 
(14) 



5 ± J (2fcr)cos(2H)], (15) 
[6t J (2/cr)cos(2fcz)], (16) 



Jo(2fcr) = — / d4>o cos(2fcr sin</>o) (17) 
2tt Jo 

is the Bessel function of the first kind. Eqs. (|13|) -(fT6 |) 
describe the system renormalization by the HF field in 
the asymptotic limit whf — > oo. Equivalently, it can 
also be seen as the lowest order of a multiple time-scale 
formalism using the expansion parameter e — ui /whf (see 
for example [H). 
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We demonstrate theoretically and experimentally the phenomenon of vibrational resonance in a 
periodic potential, using cold atoms in an optical lattice as a model system. A high-frequency (HF) 
drive, with frequency much larger than any characteristic frequency of the system, is applied by 
phase-modulating one of the lattice beams. We show that the HF drive leads to the renormalization 
of the potential. We used transport measurements as a probe of the potential renormalization. 
The very same experiments also demonstrate that transport can be controlled by the HF drive via 
potential renormalization. 
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The control of transport is a recurrent topic in physics, 
chemistry and biology. The typical scenario corresponds 
to particles diffusing on a periodic substrate, with trans- 
port controlled by the application of dc and ac external 
fields [l|, Q . The ultimate limit for the control of trans- 
port is often the impossibility of tuning the periodic po- 
tential, as it is usually the case in solid state. 

In this work we provide a proof-of-principle of how 
this limitation can be overcome, and demonstrate the- 
oretically and experimentally the control of a periodic 
potential amplitude via a strong high-frequency (HF) os- 
cillating field. The potential is renormalized, with its am- 
plitude controlled by the strength and frequency of the 
HF field. The mechanism underlying the potential renor- 
malization is the so-called vibrational resonance, intially 
introduced Q and observed 0-0] in bistable systems. 
Our experiment uses cold atoms in a dissipative opti- 
cal lattice as a model system. However, the phenomenon 
demonstrated here is very general, and is relevant to any 
classical system of particles in a periodic potential. This 
may also offer a possibility of tuning the potential in solid 
state systems, where this is usually considered impossi- 
ble. Combined with previous work which showed how 
ac fields can be used to control transport via dynamical 
symmetry breaking [l|, Q and tunnel coupling renormal- 
ization [a Q , the present work demonstrates that a com- 
plete control of transport can be achieved via ac fields. 

Our experimental work relies on the study of the trans- 
port properties of atoms in an optical lattice for different 
strengths of the applied HF field. We will demonstrate 
that by tuning the HF field it is possible to control the 
amplitude of the potential, and to make it vanish. In 
this respect, the use of cold atoms in dissipative opti- 
cal lattices is very convenient as the transpo rt prop erties 
in these systems have been studied in detail [1 CW — ll 2I | , and 
this allows us to use the transport measurements to char- 
acterize the potential. We will provide two different sets 
of measurements, as supporting evidence of the poten- 



tial renormalization. First, we will demonstrate that the 
diffusion properties, which are known to strongly depend 
on the potential depth [Io|, G3| ; can be controlled by the 
HF field in a way which corresponds to the potential 
renormalization. Second, we will show that also directed 
transport, as induced by harmonic-mixing (HM) [13| of a 
bi-harmonic drive, can be controlled by the HF field. In 
fact, anharmonicity, together with the breaking of a dy- 
namical symmetry, leads to the creation of directed cur- 
rents in harmonic mixing. Therefore, whenever the po- 
tential, renormalized by the HF field, vanishes, directed 
transport should cease EMI. 

Before discussing the experimental results, we intro- 
duce a model useful for the understanding of the po- 
tential renormalization of a dissipative optical lattice 
as produced by a HF oscillating field. We consider 
the simplest model of a dissipative optical lattice: a 
J g = 1/2 — > J e = 3/2 atom, of mass m, illuminated 
by two counterpropagating laser fields with orthogonal 
linear polarizations. This configuration generates a ID 
optical lattice |l2j. The atom in the ± ground state ex- 
periences the potential U±(z) = Uq[— 2 ± cos(2fcz)]/2, 
where z is the laser beam propagation axis, k the laser 
field wavevector and Uq the optical lattice depth [12, QjJ ■ 

We now introduce a HF oscillating force with frequency 
whf and amplitude Ahf^ 



FhfW = ^hf sin(w H Fi + ■ 



(1) 



with </>o a (mainly irrelevant) phase which describes the 
state of the oscillating force at t = 0. Of interest here is 
the high-frequency case, where the frequency of the HF 
drive is much larger then any characteristic frequency of 
the system, in the present case the vibrational frequency 
to v of the atoms at the bottom of the well. In the asymp- 
totic limit of infinite amplitude and frequency of the drive 
(whf — > oo, ^4rf — > oo), it is p ossible to show [18| that, 
consistently with Refs. [3l Il44l6| . the atomic dynamics 
corresponds to the motion in a static (i.e. without HF 



2 



field) dissipative optical lattice, with renormalized am- 
plitude U : 



U±(z) = U [-2 ± J (2fcr) cos(2fcz)]/2 



(2) 



where Jo is the Bessel function of the first kind, and 
r = Ahf/(™^hf) i s the parameter -here and thereafter 
termed the HF ratio- which controls the renormalization 
of the optical lattice. 

The above analysis shows that, in the asymptotic limit 
of infinite frequency and strength, a HF field leads to an 
effective renormalization of the potential. We now con- 
sider finite values of driving away from infinity that are 
experimentally accessible. The atomic transport in the 
optical lattice in the presence of a HF field is numerically 
studied for two different set-ups, which correspond to the 
ones used to provide the experimental evidence. 

In the first set-up, a HF force of finite amplitude and 
frequency is applied to atoms in a dissipative optical lat- 
tice. For this scheme, the effective renormalization of 
the optical potential can be detected by studying the dif- 
fusion properties of the atoms through the lattice. In 
fact, for a dissipative optical lattice of the type consid- 
ered here, it is well established [H| that there is a 
critical potential depth located at about U CT ~ 100E r , 
with E r = h 2 k 2 /(2m) the recoil energy, which sepa- 
rates two very different regimes. For potential depths 
larger than the critical one, the diffusion is normal. In- 
stead, for potential depth lower than the critical one, 
the diffusion becomes anomalous, with the exponent of 
the diffusion dependent on the potential amplitude. To 
be quantitative, we define the diffusion exponent a as 
(x 2 (t)} — (x(t)) 2 ~ t a in the limit t — > oo. According 
to this definition, a = 1 corresponds to normal diffusion, 
while a > 1 characterizes superdiffusion. In an undriven 
optical lattice, superdiffusion is encountered at shallow 
optical potentials (below the critical depth), with the ex- 
ponent a increasing for decreasing potential depth. Thus 
we will take the exponent of the diffusion as a measure of 
the potential depth. To assess the effective renormaliza- 
tion of the potential by a HF field, we numerically simu- 
lated the dynamics of the atoms in a deep optical lattice 
with a HF drive. We determined the diffusion exponent 
as a function of the HF ratio r which for infinite fre- 
quency and amplitude of the drive determines the lattice 
renormalization. Our results, reported in Fig. [TJ show 
that the diffusion exponent a can be controlled by the 
HF field, with a dependence consistent with the poten- 
tial renormalization derived in the infinite limit, Eq. ©: 
a increases whenever the potential depth is decreased, 
with the largest values of a produced by the r values 
corresponding to the zeros of the Bessel function, i.e. to 
vanishing potentials. The upper bound of a = 3 for 
Uq = corresponds to the vanishing of the friction mech- 
anism (Sisyphus cooling [l^]) associated with the optical 
lattice. Figure [1] also reports the value of Uq which cor- 
responds to the exponent a for an undriven lattice, so to 
make explicit the correspondence between a driving with 
HF ratio r and the depth of the renormalized potential. 



Finally, we notice that our results for the diffusion expo- 
nent a essentially coincides with the values derived in the 
infinite limit. We can thus conclude that the potential 
is effectively renormalised according to the dependence 
obtained in the infinite limit (see Eq. ([2])). 




FIG. 1. Left axis: numerical results, as obtained by Monte 
Carlo simulations, for the spatial diffusion exponent as a func- 
tion of the HF ratio r for an optical lattice with a depth 
Uo = 200-Er. Right axis: value of Uo which corresponds to 
the exponent a for an undriven lattice. Triangles correspond 
to results obtained in the infinite limit, for a photon scattering 
rate V' = 5lu t and T' = 10w r , were uj r is the recoil frequency. 
The filled diamonds refer to simulations with a HF field of fi- 
nite amplitude, with frequency Whf/w„ = 20 and Y 1 — 10oi r . 
The solid line is a guide for the eye for the results correspond- 
ing to the infinite limit with V' = I0uv- The dotted line is 
\Jo(2kr)\. The error bars on the numerical results correspond 
to the finite statistics of the Monte Carlo simulations. 

In the second set-up, besides the HF drive, a bi- 
harmonic force of the form 



F(t) = F L4i cos(wt) + A 2 cos(2wi + </>)] 



(3) 



is also applied to the atoms in the lattice. Here u) is 
the frequency of the drive, of the same order of mag- 
nitude or smaller than the vibrational frequency, and <p 
the relative phase between harmonics. The amplitude of 
the lattice, and its renormalization by the HF field, can 
be determined by observing the directed motion of the 
atoms through the lattice. In fact, the two harmonics of 
the drive are mixed by the nonharmonic potential, thus 
producing directed motion of the atoms through the lat- 
tice [191 ]. The average current being proportional to the 
nonharmonicity of the potential, directed transport mea- 
surements give access to the potential amplitude. More 
precisely, for weak driving the average atomic velocity is 
expected to be of the form v — v max sin(0 — <p$) , with <fid 
a dissipation-induced phase lag [19| . In our simulations 
we determined the velocity v for different values of the 
phase (j), so to derive the maximum velocity v ma x- Then 
by varying the strength of the HF drive, we were able to 
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determine v max as a function of the r-parameter, with re- 
sults as in Fig. Once again, the results produced with 
a field of large, but finite, frequency and amplitude essen- 
tially coincide with those obtained in the infinite limit. 
These results also show that current measurements can 
be used to probe the potential renormalization. When- 
ever the HF field leads to a shallower potential, as from 
Eq. @, the current is reduced, with zero current ob- 
served for those values of r leading to a vanishing poten- 
tial. 
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FIG. 2. Numerical results, as obtained by Monte Carlo sim- 
ulations, for the amplitude of the current u m ax, rescaled by 
the recoil velocity function of the HF ratio r un- 

der a biharmonic driving force of the form of Eq. © with 
Ax = A2 = 1, fb = 140hku) r and u = u v . The solid line 
corresponds to the results obtained in the infinite limit and 
the diamonds to the simulation results with ojhf/^v = 20, 
both cases with T' = 10oj r . The dotted line is |Jo(2fer)|. The 
error bars on the numerical results correspond to the finite 
statistics of the Monte Carlo simulations. 

Our experimental demonstration of potential renor- 
malization via HF field relies on the two detection 
schemes outlined above. In both set-ups, 87 Rb atoms are 
cooled and trapped in a magneto-optical trap (MOT). 
After a compression phase of 50 ms, and 8 ms of opti- 
cal molasses, the atoms are loaded into a ID dissipative 
optical lattice. The lattice is created by the interfer- 
ence of two linearly polarized and counter-propagating 
laser beams, red detuned from resonance with the D2- 
line F g = 2 — > F g = 3 atomic transition. One of the lat- 
tice beams is sent through a double pass electro-optical 
modulator (EOM), so to be able to apply a HF phase- 
modulation. In the reference frame of the lattice, such a 
phase modulation translates into a rocking force of the 
form of Eq. ([1]). Quantitatively, a phase modulation a(t) 
leads to a force F(t) — ma(t) / (2k) in the reference frame 
of the lattice In the experiments, the modulation is 
progressively turned on starting after 1 ms equilibration 
time in the optical lattice, with a turn-on ramp of 1 ms. 
Thereafter the procedure differed for the two set-ups. 

In the first experiment, we study the diffusion of the 



atoms in the optical lattice in the presence of the HF 
drive. The width of the atomic cloud is measured by 
fluorescence imaging after diffusive expansion inside the 
driven lattice. The width is measured at a fixed set of ex- 
pansion times within the interval between 1 ms and 16 ms 
from the lattice turn-on. The time range over which im- 
ages are taken is limited by the atom loss, particularly im- 
portant for the values of the HF ratio r leading to a van- 
ishing renormalized potential. Measurements of the spa- 
tial width of the atomic cloud on such a short temporal 
range do not allow us to derive an accurate value of the 
exponent of the diffusion a (2fj| . Instead, we characterize 
the diffusion by an effective diffusion coefficient D, as ob- 
tained by fitting the data with (x 2 (t)) - (x(t)) 2 = 2Dt. 
Clearly, superdiffusion leads to a large enhancement of 
the derived effective diffusion coefficient. Thus a large 
increase in the effective diffusion coefficient can be taken 
as signature of the reduction of the potential, as produced 
by the renormalization by the HF field. 
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FIG. 3. Experimental results for the effective diffusion coeffi- 
cient D as a function of the HF ratio r for different values of 
ojhf, as indicated in the figure. The data are rescaled by the 
value of the diffusion constant for an undriven lattice. The vi- 
brational frequency of the atoms at the bottom of the well, as 
determined by measuring the lattice beam power and waist, 
is u> v = (9 ± 1) • 10 5 rad/s. The solid line, with values on the 
right axis, is \ Jo(2kr)\. 

Our experimental results for the effective diffusion co- 
efficient as a function of the HF ratio r are reported in 
Fig. [3l The data clearly show that the atomic diffusion 
is significantly modified by the HF drive, with a depen- 
dence of the effective diffusion coefficient on the HF ratio 
r consistent with the potential renormalization (see e.g. 
Eq. ([2]) for the analytic expression in the limit of infinite 
frequency and amplitude). Indeed, the effective diffusion 
coefficient increases whenever the HF ratio r corresponds 
to decreasing depth of the optical lattice, with the largest 



4 



values of the diffusion constant observed in correspon- 
dence of the values of r leading to a vanishing (in the in- 
finite limit) optical lattice. This shows that the HF drive 
renormalizes the optical potential, in agreement with the 
general theory [3] and with our numerical analysis for the 
specific system. 

In the second experiment, we probe the amplitude of 
the renormalized potential by studying directed trans- 
port following harmonic mixing of two harmonics, as out- 
lined in the numerical analysis. With respect to the pre- 
vious experiment devoted to the study of the atomic dif- 
fusion, an additional bi-harmonic drive, with frequencies 
cj, 2oj and phase difference </> is introduced. This is done 
using additional acousto-optical modulators (AOMs) . In 
the reference frame of the lattice, the bi-harmonic phase 
modulation corresponds to a driving force of the form of 
Eq. ([3]). In the experiment, the HF driving is first ramped 
up, as in the previous experiment. Then the biharmonic 
drive is progressively turned on with a ramp-up time of 
4 ms. The velocity of the center-of-mass of the atomic 
cloud is derived by position measurements obtained via 
fluorescence imaging. The measurements are repeated 
for 10 different values of the phase difference cj> between 
harmonics. The data are then fitted by the expected de- 
pendence v = w max sin(0 — /fid), thus deriving a value for 
f max which can be taken as a measure of the renormal- 
ized potential depth. In fact, the mixing of harmonics 
requires an anharmonic potential, with the current gen- 
erated proportional to the anharmonicity. Our results 
for i> max as a function of the HF ratio r are presented in 
Fig. |4] These data for the directed transport amplitude 
are consistent with the renormalization of the potential 
depth by the HF drive. In fact, whenever the value of 
the HF ratio r corresponds to a reduced potential depth, 
the current decreases, with zero current observed for the 
r-values corresponding to the zeros of the Bessel func- 
tion, a signature of the vanishing optical lattice. These 
results also demonstrate a new scheme for the control of 
the transport via ac fields: the amplitude of the current 
can be controlled by a variation in the HF field and the 
direction reversed via a 7r-shift in the relative phase be- 
tween harmonics. Finally, we notice that there is a small 
deviation, both in the experiment and in the numerical 
simulations (see Fig. [5]) from the behaviour expected from 
the Bessel function at small values of r, with the data 
showing an extra peak at kr ~ 0.75. This peak could 
be explained by a superimposed resonance correspond- 
ing to the matching of the frequency of the biharmonic 
force with the oscillation frequency of the atoms at the 
bottom of the renormalized well. 



In conclusion, in this work we demonstrated experi- 



mentally the phenomenon of vibrational resonance in a 
dissipative optical lattice. The application of a HF drive, 
with frequency much larger than any characteristic fre- 
quency of the system, leads to the renormalization of the 
potential. The renormalized amplitude can be controlled 
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FIG. 4. Experimental results for the amplitude of the current 
«max, rescaled by the recoil velocity v r , of directed transport 
through the optical lattice as a function of the HF ratio r. 
In addition to the HF drive, a biharmonic force of the form 
of Eq. ([3]) is applied to the atoms, with parameters Ai = 1, 
A 2 = 2, uj = 9.42 • 10 4 rad/s, F = U2Kku r . The vibrational 
frequency of the atoms at the bottom of the well is cj v = 
(9 ± 1) • 10 J rad/s. The solid line, with values on the right 
axis, is |Jo(2fcr)|. 

by the HF drive parameters. We used transport measure- 
ments as a probe of the potential renormalization. The 
very same experiments also demonstrated that transport 
can be controlled by the HF drive via potential renor- 
malization. 

The possibility to renormalize a potential via ac fields, 
as demonstrated here, is very general, and it is applica- 
ble to any system of particles in a periodic potential. As 
such, it paves the way to the control of potentials in sys- 
tems in which they are not directly accessible, and it may 
also be applicable to solid state systems where ac drives 
can be introduced by the application of electric fields. 

Finally, our set-up can also be taken as the demon- 
stration of a sensor able to detect signals with fr equency 
exceeding any internal frequency of the sensor [14H16| . 
Here, the signal detected is the HF drive whose presence, 
although not coupling to any internal mode of the sys- 
tem, can be precisely detected due to its effect via the 
potential renormalization. 
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